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Abstract 



When a composite containing nonlinear dielectric particles suspended in 
a host medium is subjected to a sinusoidal alternating current (ac) electric 
field, the dielectric response of the composite will generally consist of ac fields 
at frequencies of higher-order harmonics. For an electrorheological (ER) solid 
under structure transformations due to external fields, we apply the Ewald- 
Kornfeld formulation to derive the local electric fields and induced dipole mo- 
ments explicitly, and then we perform the perturbation expansion method to 
extract their fundamental and third-order harmonics analj^ically. It is shown 
that the degree of anisotropy of the ER solid can affect these harmonics signif- 
icantly. Our results are well understood in the spectral representation theory. 
Thus, it seems possible to perform a real-time monitoring of the structure 
transformation by measuring the nonlinear ac responses of ER solids. 
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I. INTRODUCTION 



Electrorheological (ER) fluids [1-4] have received much attention due to the potential 
application of the rapid field-induced aggregation [2] and the large anisotropy [5], and they 
are actually a suspension in which the induced dipole moment can order the suspended 
polarizable particles into columns under the application of a strong electric field. Tao and 
Jiang examined the temporal evolution of three-dimensional structure in ER fluids by using 
a computer simulation [6] . For a wide range of the ratio of the Brownian force to the dipolar 
force, the ER fluid was found [6] to have a rapid chain formation followed by aggregation of 
chains into thick columns which has a body-centered tetragonal (bet) lattice structure, and 
the Peierls-Landau instability of single chains helps formation of thick columns. 

An ER fluid can turn into an ER sohd as the external field exceeds a critical field. For 
the ER solid, it is known that its ground state is a bet lattice. One [7] proposed that 
a structure transformation of the ER solid from the bet lattice to some other lattices can 
appear when a magnetic field is simultaneously applied perpendicular to the electric field and 
the particles have magnetic dipole moments. Then, this proposal was verified experimentally 
and a structure transformation from the bet lattice to the face-centered cubic (fee) lattice 
was observed indeed [8]. Recently, Lo and Yu have also shown that an alternative structure 
transformation from the bet structure to the fee can appear under the application of electric 
fields only [9]. 

Finite-frequency responses of nonlinear dielectric composite materials have also attracted 
much attention both in research and industrial applications during the last two decades [10]. 
In particular, when a composite containing nonlinear dielectric particles embedded in a linear 
(or nonlinear) dielectric host medium is subjected to a sinusoidal alternating current (ac) 
electric field, the dielectric response of the composite will, in general, consist of ac fields 
at frequencies of higher-order harmonics [5,11-16]. In experiments, a convenient method of 
probing the nonlinear characteristics of the composite is to measure the harmonics of the 
nonlinear polarization under the application of a sinusoidal electric field [2]. In this case, the 
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strength of the nonhnear polarization should be reflected in the magnitude of the harmonics. 
Prom the theoretical point of view, the perturbation expansion [5,14,16] and self-consistent 
methods [5,17] can be used for extracting such kind of harmonics. 

In this work, to investigate the lattice effect on the nonlinear ac responses, we shall use the 
Ewald-Kornfeld formulation [9,18-20] to derive the local electric fields and induced dipole 
moments in the ER solid which is subject to a structure transformation due to external 
fields, and then perform the perturbation expansion method to obtain their fundamental 
and third-order harmonics. To this end, it is shown that the fundamental and third-order 
harmonics are sensitive to the degree of anisotropy of ER solids. Thus, by measuring the 
nonlinear ac responses, it seems possible to perform a real-time monitoring of the structure 
transformation of ER solids. 

This paper is organized as follows. In Sec. II, we apply the Ewald-Kornfeld formula- 
tion [9,18,19] to derive the local electric field and induced dipole moment in ER sohds, and 
then perform the perturbation expansion method to extract their fundamental and third- 
order harmonics. In Sec. Ill, we numerically investigate these harmonics under different 
conditions, which is followed by a discussion and conclusion in Sec. IV. 

II. FORMALISM 

A. Local electric field and induced dipole moment: Ewald-Kornfeld formulation 

Let us start by considering the ground state of an ER solid, namely a bet (body-centered 
tetragonal) lattice, which can be regarded as a tetragonal lattice, plus a basis of two particles 
each of which is fixed with a point dipole at its center. One of the two particles is located at 
a corner and the other one at the body center of the tetragonal unit cell. Its lattice constants 
are denoted by a{— b) — iq"^^'^ and c — qi along the x{y) and z axes, respectively. In this 
case, the uniaxial anisotropic axis is directed along the z axis. As q varies, the volume of the 
unit cell keeps the same, Vc = i^. Consequently, the degree of anisotropy of the tetragonal 
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lattice is measured by how q is deviated from unity. 

— * 

In case of an x— directed external electric field £'0, the dipole moments p — px are 
perpendicular to the uniaxial anisotropic axis. Then, the local field E (e.g., E = E^x, 
Ez — 0) at the lattice point ^ = has the following form [9,18-20] 

E.-pj:j:[-j,{R,)+xyj,{R,)] - ^ E n(G)§exp(^) + f^. (1) 
In this equation, 71 and 72 are two coefficients, given by 

3erfc(?7r) Arf 677 
72 r = + {-1=^ + -J=-^) eM-V r ), 

where erfc(?7r) is the complementary error function, and r] an adjustable parameter making 
the summation converge rapidly. In Eq. (1), R and G denote the lattice vector and the 
reciprocal lattice vector, respectively, 

R = e{q-y^lx + q-^/'^my + qnz), 
G — —^{q^^'^ux + q^^'^vy + q~^wz), 

where /, m, n, u, v, w are integers. In addition, Xj and Rj of Eq. (1) are given by, respectively. 



Xj = I — — , Rj = 1^ — — {ax + ay -\- ci)|. 



and the structure factor n(G) = 1 + exp[i(M + v + w)/t^]. 

So far, based on a self- consistent method, we apply the result of the local field to evaluate 
the effective polarizability cte// of the dipole lattice, 

a , , 

where a stands for the polarizability of an isolated dipole, and — E^Vc/p the local field 
factor which was measured in the computer simulations by Martin et al. [21,22]. Let us 

use 7^ and 7^; (= 7^,) to denote the local-field factors parallel and perpendicular to the 
uniaxial anisotropic axis, respectively. Accordingly, we have 7^ = 7^ = 47r/3 for the bcc 



lattice {q — 1). In what follows, we set 7' = 3^ /An. It is worth noting that 7' is a function 
of a single variable, namely degree of anisotropy q. Also, 7^ and 7^ satisfy the sum rule 
7^ + 27^ = 3 [23], and 7^ = 7^ = 1 [at g = 1 (bcc)] just represents the isotropic limit. 

Based on Eq. (2), it is straightforward to derive the dipole factor (also called Clausius- 
Mossotti factor) for a specific particle b'^ and b'^ along the x (or y) and z axes, respectively 

where p stands for the volume fraction of the particles, and b the known dipole factor for an 
isolated particle, b = (ci — e2)/(ei + 2e2). Here ei and 62 denote the dielectric constants of the 
particle and host fluid, respectively. The two parameters, b'^ and b'^, have been represented 
as by using the joint subscript x{z). The similar notation are also used for 7^.(^). In what 
follows, for convenience, will be used to stand for X^"^^ and X^^^ which represent the 

transverse [i.e., the external field is parallel to the x (or y) axis] and longitudinal [namely, 
the external field is parallel to the z axis] field cases of the quantity X, respectively. 

So far, the local-field effect arising from all the other particles in the lattice has been 
exphcitly included in Eq. (3) by using the Ewald-Kornfeld formulation. Next, let us see the 
particle as the one embedded in an effective medium with an effective dielectric constant 
which can be determined by 

It is worth remarking that this e^^^^^^^ has included the detailed structural information of 
the lattice, as expected. 

— * 

As the external field Eo is along the x or z axis, the volume average of the local electric 
field is given by 

where (■ ■ ■) stands for the volume average of • • •. In view of Eqs. (3) and (4), e^"^-*'^^-*! of 
Eq. (5) is given by 



^'P7i(.) - 26 - 1 " 

Then, it takes one step forward to express the corresponding induced dipole moment 
inside the particle 

^TML, ^ epm^rX^^^Eo. (7) 

So far, we have derived the local electric fields [Eq. (5)] and induced dipole moments 
[Eq. (7)], by taking into account the lattice effect in a rigorous manner. In what follows, 
we shall use pi (or pi), 7', Ei (or £"1) and to denote both the transverse and longitudinal 
field cases, if there are no special instructions. 



B. Nonlinear polarization and its higher-order harmonics 

For studying the nonlinear ac responses, let us consider the particle with a cubic nonlin- 
earity like 

ei = ei + xi£^i ~ei + Xi(£^i'), (8) 

where Xi represents the nonlinear susceptibility of the particles embedded in a linear host 
fluid. Throughout the paper, we focus on the weak nonlinearity only, namely xii^i^) ^ ^i) 
and the low concentration limit, and we use "~" to indicate the quantities having a nonlinear 
characteristic. 

If a sinusoidal electric field 

Eo{t) ^ Eo sin ujt (9) 

is applied, the local electric field (£'1) and the induced dipole moments pi will depend on 
time sinusoidally, too. 

Owing to the inversion symmetry of the system, the local electric field is a superposition 
of odd-order harmonics such that [5,14] 
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(^i) = sm{ut) + E^^ sm{?,ut) + • • • . (10) 

Accordingly, the induced dipole moment contains harmonics as [5,14] 

Pi = Poj sm{ijjt) + sin(3a;i) H . (11) 

For the purpose of normahzation, we shall use po = ^2rlbEo which is independent of the 
nonlinear characteristic as well as the degree of anisotropy 7'. 

C. Analytic solutions 

In what follows, we apply the perturbation expansion method to extract the harmonics 
of the local electric fields {E^ and E^oj) and the induced dipole moments (p^ and ps^^). 

Let us denote the volume average of the linear local electric field (namely, the local 
electric field at which the nonlinear characteristics of the particles disappears) [Eq. (5)] as 

{E^) = F{e,,e2,p,^',Eo). (12) 

In view of the existence of nonlinearity, we have 

(^i) = F(ei,e2,p,7',^o). (13) 

Next, expand the local field (£^1) into a Taylor expansion, taking Xi(-E'i) as the pertur- 
bative quantity. As a result, we obtain [5] 

(^1) = F(ei, £2, P, 7', ^0) + -^/i^u £2, P, 7', Eo)U,=e,Xi{E'i) + ■■■■ (14) 

Keeping the lowest orders of xi^oit)'^ and xi{Ef), because of the time-dependence of the 
external electric field [Eq. (9)], we obtain the local electric fields and induced dipole moments 
in terms of the harmonics, E^^, £'3^^, p^^, and psuj [note the higher-order harmonics (namely, 
fifth-, seventh-order, et al.) have been neglected]. Here the harmonics of the local electric 
fields and induced dipole moments are respectively given by 
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where 



with 



xl^'E,^ - ^(x}/'^o)^ (16) 



4W. 



4Vr2^e2(ei - €2) 

Po = . U-f ^ 7^ 2^ <9iv Xi ^0) , (18) 

462^2^ ei - e2 7'e2P+ ei 3 - Yp 2 ^ 



W^i = -7'eiP + e2(3 + 7V), 

1^2 = ei - 7'eip + 62(2 + 7'p). 
4eil^2(ei - 62) 



7'e2P + ei(3 - 7'p) ' 
3Wi 



- [7'eip(2 + yp)Vc/, + C/2]2 E Q-(x}/^£^o)^ 



gi=4ye?e2p(3 + 7'p-7V), 
g2=7'eiP(3-47'p + 7V), 
<53 = 7'e2P(6 + 57'p + 7V), 
Q4 = -27'eie^p(6 + 77'p + 27'V'), 

= 3e2e2(-9 - 37'p + 37'^' + 27'V'), 



III. NUMERICAL RESULTS 

We are now in a position to do some numerical calculations to discuss the effect of 
the degree of anisotropy q on the harmonics of the local electric fields and induced dipole 
moments. Without loss of generality, take ei = 30eo, £2 = 2.25eo (dielectric constant of 
silicone oil), and p = 0.1 for numerical calculations, where eo denotes the dielectric constant 
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of free space. In fact, here p — 0.1 corresponds to the dilute hmit. Nevertheless, it would 
be reahzable as the particles can have a (solid) hard core and a relatively soft coating (by 
long polymer chains, etc.) to avoid aggregation. Thus the model of a soft sphere with 
an embedded point dipole can be used. In this case, the many-body (local field) effects 
are important, while the multipolar effects can be neglected. This is exactly part of the 
emphasis in the present work. On the other hand, for a higher p, similar results can be 
achieved indeed, as expected. 

In Fig. 1, we investigate the dependence of the local- field factor [4] on the degree of 
anisotropy q. Increasing q causes 7^ (or 7^) to increase (or decrease). Also, a plateau is 
shown at 7^ = 7^ = 1, which actually includes the transformations ranging from the bcc 
(g — 1) lattice to the fee {q — 2^/^). Accordingly, similar plateau occurs at all other figures 
(see Figs. 2 and 3). 

Figure 2 displays the fundamental and third-order harmonics of the local electric field 
and induced dipole moment, as a function of the degree of anisotropy q for different Xi'^Eq 
for the transverse field case. It is found that increasing q causes both the fundamental and 
third-order harmonics of the local electric field to decrease. Accordingly, the fundamental 
harmonics of the induced dipole moment is caused to decreasing as well. However, the 
third-order harmonics of induced dipole moment is caused to increase as q increases. Note 
the three curves in Fig. 1(c) are overlapped. 

Similarly, we also investigate the harmonics of the local electric field and induced dipole 

moment in Fig. 3, but for the longitudinal field case. In contrast to the transverse field case 

(Fig. 2), the exactly opposite effects have been shown for the longitudinal field case (Fig. 3), 

due to the opposite dependence of q on 7^ and 7^, see Fig. 1. Again, the three curves in 

Fig. 1(c) are overlapped. 

Both Figures 2 and 3 show that the harmonics are strongly dependent on the nonlinear 

1/2 

response of the suspended particles. Moreover, increasing the nonlinear response X\ -E'o 
leads to an increase in the harmonics. This is in agreement with the results of our recent 
work [5,24]. In addition, we also studied the effect of volume fraction p, and found that 
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increasing the volume fraction can enhance the harmonics too (no figures shown here). 

To understand the results, we resort to the spectral representation approach [25]. Let's 
start by denoting s = (1 — ei/e2)~^, and then the local electric field given by Eq. (5) admits 

(El) = Eo + -^Eo, (19) 

S — Si 

where the residue Fi = 1/3 and the pole si = (1 — pj')/3. As p — > 0, si tends to 1/3, which 
produces the known value for an isolated particle. The substitution of the present model 
parameters yields s — —0.081, and hence we have 

<^-> = ^o+ -0.081- If- Q.iy)/3 ^°- 

Based on this equation, it is apparent to see that increasing 7^ or 7^ causes the electric field 
to decrease. In addition, we have shown that increasing q causes 7^ (or 7^) to increase (or 
decrease), see Fig. 1. Also, a plateau was shown at 7^ = 7^ = 1, which actually includes the 
transformations ranging from the bcc {q — 1) lattice to the fee (g — 2^/^). Because of the 
dependence of the harmonics on 7', similar plateau have already been shown in Figs. 2 and 3, 
too. As a matter of fact, the effective hnear dielectric constant mainly contributes to the 
magnitude of the fundamental harmonics of the induced dipole moment or local electric field. 
In the meantime, because of the existence of weak nonlinearity, the nonlinear polarization 
just has a perturbation effect, which can be neglected when compared to the linear part. 
Since the spectral representation approach is a different way of formulating the linear local 
electric field, to some extent, it can be used to explain the fundamental harmonics of the 
local electric field and hence those of the induced dipole moment which are displayed in 
Fig. 2(a,c) and Fig. 3(a,c). Based on Eq. (20) (namely, an expression obtained from the 
spectral representation approach), it is convenient for us to understand that a decrease 
(or increase) of q should lead to an increase of the local electric field for transverse (or 
longitudinal) field cases. In this connection, for transverse (or longitudinal) field cases, the 
fundamental harmonics of the local electric field and hence corresponding induced dipole 
moment should increase as q decreases (or increases). This has also been displayed in 

10 



Fig. 2(a,c) (or Fig. 3(a,c)). The explanation on the third-order harmonics of the local 
electric field and induced dipole moment seems to be more complicated. Nevertheless, since 
the strength of the nonlinear polarization of nonlinear materials should be reflected in the 
magnitude of the high-order harmonics of local electric fields and induced dipole moments, 
the effective nonlinear part in the effective nonlinear dielectric constant should be expected 
to determine the magnitude of the third-order harmonics. It appears that the third-order 
harmonics of the local electric field show a behavior similar to its fundamental harmonics. 
However, inverse behavior occurs to the third-order harmonics of the induced dipole moment. 
This difference would possibly be due to the nonlinearity fluctuation which might arise from 
the long-range interaction. This interaction existing in the system is explicitly evaluated in 
the current paper, by using the Ewald technique. 



IV. DISCUSSION AND CONCLUSION 

Here some comments are in order. In this work, we have used the Ewald-Kornfeld 
formulation to derive the local electric fields and induced dipole moments in a rigorous 
manner, in an attempt to take into account the lattice (local-field) effect on the nonlinear 
ac responses of ER solids under structure transformations. 

In fact, we presented an (approximate) effective medium theory (EMT) for considering 
the local-field effect on the electrorotation and dielectric dispersion spectra of colloidal par- 
ticles or biological cells in the previous paper [26]. However, we were unable to study the 
detailed structural information via the EMT. In the present paper, we have considered the 
local-field effect in a rigorous manner, based on the Ewald-Kornfeld formulation which is 
quite different from the EMT. 

We have considered the fundamental and third-order harmonics. As a matter of fact, 
it should be straightforward to investigate the higher-order (e.g., fifth-order, seventh-order, 
etc.) harmonics by keeping more terms in Eq. (14). Due to the virtue of symmetry, the 
odd-order harmonics appear only. In fact, the even-order harmonics can also be induced 
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to occur due to the coupling between the ac and dc external fields [16]. Since the second- 
order harmonics can often be of several orders of magnitude larger than the third-order, 
it is of value to discuss the lattice effect on the nonlinear ac responses by applying an ac 
field coupling with a dc field, in order to investigate the second-order harmonics. Finally, 
this work can be extended to poly disperse ER solids [27] in which the permittivities of 
particles can possess a distribution, or to a system [28] where anisotropics of material and 
geometry are caused to exist due to fabrication methods or field effects. On the other hand, 
the nonlinear response of composites to ac fields will become more complicated when the ac 
field contains different basic frequencies or different magnitudes of external fields at different 
basic frequencies [29] . 

Our former work [5] covers a wider range of the degree of anisotropy q. However, the 
explicit form of the local field-g relation was lacking. By invoking a lattice structure (which 
exists in ER sohds), this can be achieved, as evident in this work. In detail, in Ref. [5] 
we used an effective medium theory to derive effective (anisotropic) dielectric constants, 
based on the Maxwell-Garnett approximation. In contrast, in the current paper, the Ewald 
technique is used to derive the effective dielectric constants, by including the long-range 
interaction explicitly (which was, however, lacking in Ref. [5]). In this sense, for treating 
ER solids with lattice structures, the present derivation of the effective dielectric constants 
should be expected to be more accurate. Interestingly, some result of the current paper was 
also predicted in Ref. [5] . For instance, as q decreases, for longitudinal (or transverse) field 
cases 7^ (or 7^) is caused to increase (or decrease), thus yielding an increasing (or decreasing) 
third-order harmonics of the induced dipole moment, see Fig. 3(d) (or Fig. 2(d)). Here 7^ 
(or 7^) corresponds to /3l (or /3t) of Ref. [5]. It is worth to mentioning that a different 
normalization factor (i.e., po) was used in Ref. [5]. 

To sum up, based on the Ewald- Kornf eld formulation, we have investigated the nonlinear 
as responses of an ER solid with nonlinear spherical particles embedded in a linear host fluid, 
and found that the fundamental and third-order harmonic ac responses are sensitive to the 
degree of anisotropy within the ER solid. Our results have been well understood through 
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the spectral representation theory. Thus, by measuring the nonhnear ac responses of ER 
sohds, it seems possible to perform a real-time monitoring of structure transformations. 

ACKNOWLEDGMENTS 

This work was supported by the Alexander von Humboldt foundation in Germany. The 
author acknowledges Professor K. W. Yu's very fruitful collaboration and comments. 



13 



REFERENCES 

[1] Halsey, T. C. Science 1992, 258, 761. 

[2] Klingenberg, D. J. MRS Bull 1998, 23, 30. 

[3] Wen, W.; Huang, X.; Yang, S.; Lu, K.; Sheng, P. Nat. Mater. 20 03, 2 , 727. 

[4] Huang, J. P. Phys. Rev. E 20 04, 70, 042501; ihid. Chem. Phys. Lett. 2004, 390, 380. 

[5] Huang, J. P.; Wan, J. T. K.; Lo, C. K.; Yu, K. W. Phys. Rev. E 2001, 64, 061505(R). 

[6] Tao R.; Jiang, Q. Phys. Rev. Lett. 1994, 73, 205. 

[7] Tao, R.; Jiang, Q. Phys. Rev. E 1998, 57, 5761. 

[8] Wen, W.; Wang, N.; Ma, H.; Lin, Z.; Tarn, W. Y.; Chan, C. T.; Sheng, P. Phys. Rev. 
Lett. 1999, 82, 4248. 

[9] Lo, C. K.; Yu, K. W. Phys. Rev. E 2001, 64, 031501. 

[10] Bergman D. J.; Stroud, D. Solid State Phys. 1992, 46, 147. 

[11] Levy, O.; Bergman, D. J.; Stroud, D. Phys. Rev. E 1995, 52, 3184. 

[12] Hui, P. M.; Stroud, D. J. Appl. Phys. 1997, 82, 4740. 

[13] Hui, P. M.; Cheung, P. C; Stroud, D. J. Appl. Phys. 1998, 84, 3451. 

[14] Gu, G. Q.; Hui, P. M.; Yu, K. W. Phystca B 2000, 279, 62. 

[15] Huang, J. P.; Gao, L.; Yu, K. W. J. Appl. Phys. 2003, 93, 2871. 

[16] Wei, E. B., Song, J. B.; Gu, G. Q. J. Appl. Phys. 2004, 95, 1377. 

[17] Wan, J. T. K.; Gu, G. Q.; Yu, K. W. Phys. Rev. E2001, 63, 052501. 

[18] Ewald, P. P. Ann. Phys. (Leipzig) 1921, 64, 253. 

[19] Kornfeld, H. Z. Phys. 1924, 22, 27. 

14 



[20] Lo, C. K.; Wan, J. T. K.; Yu, K. W. J. Phys.: Condens. Matter 2001, 13, 1315. 

[21] Martin, J. E.; Anderson, R. A.; Tigges, C. P. J. Chem. Phys. 1998, 108, 3765. 

[22] Martin, J. E.; Anderson, R. A.; Tigges, C. R J. Chem. Phys. 1998, 108, 7887. 

[23] Landau, L. D.; Lifshitz, E. M.; Pitaevskii, L. P. Electrodynamics of Continuous Media, 
2nd ed.; Pergamon Press: New York, 1984. 

[24] Huang, J. P.; Gao, L.; Yu, K. W.; Gu, G. Q. Phys. Rev. E 2004, 69, 036605. 

[25] Bergman, D. J. Phys. Rep. 1978, 43, 377. 

[26] Gao, L.; Huang, J. P.; Yu, K. W. Phys. Rev. E 2003, 67, 021910. 
[27] Sun, H.; Yu, K. W. Phys. Rev. E 2003, 67, 011506. 
[28] Huang, J. P. J. Phys. Chem. B 2004, 108, 13901. 

[29] Wei, E. B.; Yang, Z. D.; Gu, G. Q. J. Phys. D: Appl. Phys. 2004, 37, 107. 



15 



FIGURES 

FIG. 1. Local field factor 7' versus the degree of anisotropy q. The bet (ground state), bcc and 
fee lattices which are respectively related io q = 0.87358, q = 1.0 and q = 2^/^ are also shown 
(dot-dashed lines). 

FIG. 2. The fundamental and third-order harmonics of the local electric field and induced 
dipole moment, as a function of q for Xi'^Eq = O.Te^^ (solid lines), Xi'^Eq = CSey^ (dotted 

1 /2 1 /2 

lines), and Xi ^0 = 0.9eg (dashed lines), for the transverse field ease. Note the three curves in 
(c) are overlapped. Parameters: ei = 30eo, £2 = 2.25eo, and p = 0.1. The bet, bcc and fee lattices 
are shown as well (dot-dashed lines). 



FIG. 3. Same as Fig. 2, but for the longitudinal field ease. 
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